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We study inflection point inflation using Singularity Theory, which relates degenerate critical 
points of functions to their local behavior. This approach illuminates universal features of small-field 
models and gives analytic control over parametrized families of scalar potentials near inflationary 
solutions. The behavior of the scalar potential is tied to the number of physical input parameters, 
which determines a set of universality classes. Within these classes, we obtain universal scaling 
relations for density perturbations and the scale of inflation. In specific models, we show that the 
scale of supersymmetry breaking also possesses scaling behavior. We illustrate this general structure 
with a specific example: the Racetrack Inflation model in type IIB string theory, with the inflaton 
being the real part of the Kahler modulus, and the input parameters being flux dependent quantities 
that appear in the 4D N = 1 superpotential. 



I. INTRODUCTION 

A period of inflation in the early universe, dominated by vacuum energy, solves the flatness and horizon problems of 
big bang cosmology, in addition to diluting unwanted relics |30[ I37[ 138] . The vacuum energy is converted to radiation 
after inflation ends, which acts as the source of entropy for the universe. Most importantly, quantum de Sitter 
fluctuations of the inflaton leave an almost scale-invariant spectrum of perturbations on space-time, which leads to 
the formation of structure. 

The minimum requirement for a functioning model of inflation is a scalar potential which generates sufficient e- 
foldings and is compatible with the observed values of the power spectrum Tn{k) and the spectral index n s (k). A 
successful model should also satisfy the deeper need for a consistent embedding within particle physics and string 
theory. This has been an arena of much activity; for recent reviews, we refer to |14[ 139] . 

A sensible embedding within particle physics would ostensibly address the conceptual problems associated with 
trans-Planckian VEVs of the inflaton, which are hard to understand within the framework of effective field theory. 
Notable progress in this direction has been achieved [40l [41] . 

Another interesting and much-studied route has been the class of small-field models, especially models where 
inflation occurs near an inflection point of the inflaton potential. These offer a simple way to evade trans-Planckian 
VEVs — the extreme flatness of the potential near the inflection point guarantees that enough e-foldings are generated 
while the field remains in the vicinity, provided the initial conditions are satisfied. Moreover, the scale of inflation can 
be low in these models, enabling a more direct connection with the scale of supersymmetry breaking. Inflection point 
models are common: D-brane inflation [SJ [T31 HH 122] ; Kahler modulus inflation [T^l [TTJ [23] , MSSM inflation with soft 
terms [9], and without soft terms [7] are all examples. 

From a purely mathematical perspective, inflationary models that employ inflection points reduce to the study of 
functions whose first and second derivatives nearly vanish. These inflection or degenerate critical points determine 
the qualitative behaviour of their functions. 

In this paper, we study inflection point models using degenerate critical points. The appropriate language is 
Singularity Theory which relates degenerate critical points of functions to their behavior. Since the scalar 
potential governing inflation is a real function, we apply techniques from Catastrophe Theory. Catastrophe Theory 
studies critical points of real functions and therefore naturally describes inflection point inflation. 

This mathematical language is useful in various ways. First, all models have a set of physical input parameters 
which determine the inflaton potential. Catastrophe Theory classifies the kinds of inflection points one can obtain 
by the number of control parameters. Second, for a given number of control parameters, Catastrophe Theory gives 
analytic control over the parameter space and therefore the regions where inflection points may be obtained. This 
is particularly useful for complicated potentials with many input parameters. Third, Catastrophe Theory grants 
universal scaling relations for physical observables such as energy scales and density perturbations. 

These methods are applicable to virtually all models of inflection point inflation (with caveats spelled out later) . In 
this paper we illustrate the general structure with a specific example: the well-studied Racetrack Inflation model in 
type IIB string theory [15]. The inflaton is the real part of the Kahler modulus, and the control parameters are flux 
dependent quantities that appear in the 4D Af = 1 superpotential. There is a single variable and three independent 
control parameters, corresponding to an singularity or "Swallowtail Catastrophe" . Parameter ranges which give 
acceptable inflation are obtained. 

The Hubble parameter during inflation is found to scale with the separation in field space of the point of inflation 
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and a metastable vacuum with nearly vanishing cosmological constant. The vanishingly small vacuum energy was a 
numerical challenge for previous work [8], but is now readily found and conceptually straight forward. 

For an appropriate choice of parameters, the vacuum may be widely separated from the point of inflection, leading 
to a separation of the scales of inflation and SUSY breaking [35]. Catastrophe Theory describes precisely how these 
scales separate. 

The plan of the paper is as follows. In Section [TTJ we introduce the machinery of Catastrophe Theory in a way 
that can be directly applied to inflection point inflationary models. In Section |III[ we work out a simple toy model 
of inflation with two control parameters, corresponding to the "Cusp Catastrophe" . In Section |IV| we work out a 
more realistic example of inflation with three control parameters, corresponding to the "Swallowtail Catastrophe" . 
We describe in detail the example of racetrack inflation, with particular focus on the scales of inflation and SUSY 
breaking. We conclude with a discussion. 



II. MACHINERY 



A. Morse Theory and Inflation 



The study of inflation builds on the study of the scalar potential, which determines the stability of the physical 
system. Stability hinges on the existence of real critical points of the potential — local extrema. Morse theory 
formalizes the familiar notion of stability. 

The local behavior of a function at a generic point is determined by its derivative there. Functions are typically 
increasing or decreasing. A point where the first derivative vanishes is called critical. The local behavior at a 
critical point is determined by a nonvanishing second derivative. Stability of the system can be inferred from the 
the Hessian Matrix. If all the eigenvalues are positive, the system is stable. Negative eigenvalues correspond to 
unstable directions. In scalar field theory, these would be massive and tachyonic fields, respectively. Even better, 
small perturbations of such a potential don't change its stability properties. The function is said to be structurally 
stable against perturbations. 

When an eigenvalue of the Hessian vanishes, the function is no longer described by Morse theory. The system is 
marginally stable at the critical point, but structural stability is lost. For example, x 3 has a vanishing first and second 
derivative at the origin, but a small perturbation 

x 3 + ax, (II. 1) 

is classically stable for any a < (local extrema), but unstable for any a > (imaginary critical points). As 
a approaches zero, the critical points degenerate. This non-Morse behavior of changing structural stability is the 
subject of Catastrophe Theory. 

Notice that the structure of the function in (II.l ) changes as a passed through zero. Generally, functions will depend 
on many control parameters, like a. We therefore speak of families of functions with a control parameter space, C. In 
the context of effective field theory, this is the space of couplings. As familiar from Statistical Mechanics, couplings 
can be tuned to yield Critical Phenomena, where scale invariant thermal (or quantum) fluctuations dominate the 
physical system. This occurs when critical points of the function describing the system degenerate. For inflation, 
these fluctuations manifest themselves as the nearly scale invariant spectrum of density perturbations. 

Now, a sufficient condition for inflation is that the slow roll parameters nearly vanish. Small-field inflation imposes 
a stricter version of this condition, namely, 
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(II.2) 



Thus, small-field inflation occurs near a degenerate critical point. 



B. Catastrophe Theory and Inflation 



Catastrophe Theory explains how the stability structure of Morse functions change. Such change occurs when 
critical points degenerate, which is also a condition for Inflation. The observation, illustrated by (II.l), results in 
a domain structure for C. This fact will yield important relationships between this domain structure, inflation and 
the physics in which inflationary models are embedded. This section is devoted to generalities. We will discuss two 



explicit examples of how such domain structures form in sections III and IV 
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For clarity, we begin with a single field. Let V(x) be the potential associated with a scalar field x. V is a real 
function which also depends on a number of couplings. Physically, we are often interested in the whole space of 
possible couplings, C(V). 

As discussed above, the critical points ft determine the structure of a function. To illuminate this fact, one may 
write the derivative of V as 

V'(x)=v(x)Yl(x-0i), (11.3) 

i 

where v(x) never vanishes. Clearly, the critical points ft will change with the control parameters in C(V). Since V is 
a real function, the ft may appear in complex conjugate pairs. Real roots correspond to local extrema. As discussed 
above, the stability structure of V changes (minima disappear) as ft pass between real and complex values. This 
leads to the observation of singular importance in Catastrophe Theory: such transitions necessarily occur when the 
critical points are identical — they degenerate. 

A generic point in C(V) corresponds to distinct ft, and thus a Morse function. For an n-dimensional parameter 
space, however, there exist subsets of couplings where the ft degenerate. For example, if two of the critical points 
take the same value a, 

ft = ft = a, 

this corresponds to a single constraint between the couplings in C(V). This single constraint corresponds to a 
hypersurface in C(V) of codimension one — a domain wall. Roots transition between real and complex values only 
as the couplings jump across the wall. These walls divide the space of control parameters into domains where V has 
different stability structures. 

For a generic point on the domain wall, V has a doubly degenerate critical point, a and may be expressed locally 
by its Taylor series, 

V(x) « V(a) + ^(x- af + 0[(x - a) 4 }. (II.4) 
The linear and quadratic terms vanish; V is non-Morse. More importantly, the slow roll conditions for small-field 



inflation (II. 2 1 are satisfied in, and only in, the vicinity of a domain wall. 



Higher degeneracy among the critical points can be described similarly. An m-fold degenerate critical point, 

ft = • • • = ft„ = 

corresponds to a (n — m)-dimension hypersurface in C(V). This process terminates at m = n, when the maximum 
possible critical points degenerate. This corresponds to a point C(V), c . The function, V(x), associated with the 
couplings at c is called the germ of the Catastrophe. This "domain complex" will have important implications for 
the effective field theory. 

There is a technical point worth mentioning. The mathematical space of control parameters, say C, is n dimensional, 
but C may have lower dimension. This just means that the physical parameter space is less than the most general 
allowed by mathematics. This is common, particularly given some imposed symmetry. For simplicity we shall ignore 
these cases and assumepQ C = C, and discuss explicit exclusions of parameter space where appropriate. 

Inflation occurs in the vicinity of a domain wall. Deviations of the couplings away from some point on the domain 
wall can be encoded in a nonvanishing Aj.. This parametrizes the "unfolding" of the degenerate critical point. This 



means that (II. 4 1 becomes, 

V{x) « V(a) + Xi(x - a) + ^-{x - a) 3 + 0{(x - a) 4 ]. (II.5) 
Here a is still the point of inflection: V"(a) = 0. Notice that the slow roll parameter e = ^Xf. In the slow roll 



approximation, one may compute the number of e-foldings associated with (II. 5 1 
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which, to a very good approximation, is 
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To ensure the proper numbers of e-foldings, Ai must be quite small. This leads us to our first fundamental result: 
Ai is essentially model independent and is degenerate with the number of e-foldings [2]- This is important, as A3 
deviates only infinitesimally from the domain wall value where Ai = 0. As we will demonstrate with an explicit 
example in section [TV] this distinction is physically unimportant. Thus, A3 depends solely on a subset of couplings 
which parametrize the domain wall. As we shall also see, A3 can be computed explicitly once the type of Catastrophe 
is chosen. Taken together, these facts yield our second fundamental result: all physical observables connected with 
the inflation sector are determined by the domain wall parameters, and may therefore be classified by Catastrophes. 
We now turn to the classification of these Catastrophes. 



C. ADE Classification of Inflaton Potentials 

We now present Arnold's classification of Catastrophes relevant for the scalar potential. It relies one two numbers: 
dimensionality of C(V) and number of possible vanishing eigenvalues of Hessian. There are also technically important 
minus signs. 

The Splitting Lemma is an important result for this classification. It proves the existence of a "good" coordinate 
patch on the target space E where Morse and non-Morse behavior can split into different variables. Mathematically, 
this removes any constraint on the dimensionality of E and allows us to focus only on the number of vanishing 
eigenvalues of the Hessian matrix. Physically, it means that we can always find a field redefinition which splits all 
scalar fields into an "inflaton" (or the relevant multifield configuration) and "spectators". Of course, one must always 
ensure that spectators are stabilized during inflation. 

Arnold's work involved complex functions. Since we are interested in the scalar potential, we extract what is 
relevant for real functions. Arnold demonstrated that the Catastrophe domain structures are linked to symmetries of 
the space of Taylor polynomials (Jets) at a critical point of V. Specifically, this space at a degenerate critical point 
is invariant under the action of some maximal Lie Group Q . The larger the symmetry group, the more complicated 
the behavior, and the higher order the Taylor polynomial. For explicit details of the group action see the appendix. 

Simply put, if there is only a single vanishing eigenvalue of the Hessian, Q must be the group A kl where rank k is 
the dimension of the control parameter space. If there may be two vanishing eigenvalues, the group may also be D k 
or E k . Technical issues arise if there may be more than two vanishing eigenvalues. Given such a case for dimC > 6, 
for instance, the critical point is associated with a one parameter family of germs. Q may no longer be simple, but 
may pick up a U(l) factor, for instance. While elegant, this is far beyond the scope of our work. For proofs and a 
more detailed discussion, see Arnold [TTj . 

These considerations are deep, but the results are simple. For all relevant applications, the scalar potential V is 
mapped into a finite polynomial which depends on the Catastrophe Structure of V. Usually, this map is just the 
relevant Taylor polynomial. 

In what follows, x, y are fields parametrizing the some target space E and other parameters belong to C. 

fc-i 

A ±k : (±) fe .T fc+1 + a mx m , 

m— 1 

fc-3 

D ±k : (±) k xy 2 ± x 2 ^ 1 + ^ a m x m + ay + c 2 y 2 , 

m=l 

E±q : ±(x 4 + y 3 ) + ax 2 y + bx 2 + cxy + dx + fy, 

4 

E 7 : y 3 + yx 4 + ^ + by + cxy, 

m—l 

3 3 

E s : x 5 +y 3 +yY, + CmxW ' '• 

rn—O m—l 

In this representation, the origin of C corresponds to c . 

Applying this classification to a smooth potential is straight forward. If only a single eigenvalue of the Hessian may 
vanish, then the Catastrophe depends on the number of couplings. For n couplings exclusive of a constant term, we 
have A n+ 2 behavior. When two eigenvalues may vanish, D k behavior is typically associated with dimC(V) — k. Of 
course for k = 5, 6, 7, one must be wary of exceptional germs. In this work we focus primarily on the single field case, 
leaving the more complicated multifield parameter spaces for future work. 



(II.8) 

(II.9) 
(11.10) 
(11.11) 

(11.12) 
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III. Ay. THE CUSP CATASTROPHE 



A. A3 Domain Structure 



A potential with A 3 behavior is simplest workable model of inflation. A triply degenerate critical point can unfold 
to leave a single doubly degenerate critical point for inflation. The remaining root is a local minimum - the vacuum. 
This class of potential exhibits the majority of relevant behaviors. In this section we investigate we build the domain 
structure of A3 and examine the physical details. 

The most general, local form for a potential with with A 3 behavior is 

V(x) = ^x 4 + ^ax 2 + bx + c . (III.l) 

For stability, we require A > 0. Demanding V support both inflation at a degenerate critical point as well as a stable 
vacuum elsewhere requires all three critical points of V to be real, and two to be degenerate. Inflation occurs at the 
latter, x — a, and ends with the field rolling into the final vacuum at x = j3. There the inflaton decays and reheating 
begins. In short, 

^ = \(x-a) 2 (x-0) . (111.2) 
ax 



We now illustrate some features of this simple model. Taking a derivative of Eq. III.l and matching coefficients in 
Eq. ITTPl we find 

2a + /3 = 0, (III.3) 

and the relations 

a 2 + 2a/3 = a/A, -a 2 /? = b/X . (III.4) 

Eliminating /3, one finds, 



a/A = -3a 2 , b/X = 2c? . (III.5) 

Eliminating a, we have the constraint 

(£)'+(*)'-'■ 

Eq. |IIL6] is the condition that V must satisfy in order to have an inflection point. In other words, Eq. |ffl.5] defines a 
domain wall in parameter space. The name of the Cusp Catastrophe originates from the shape of the domain wall, 
as seen in Figure [I] (with A = 1). 

It is instructive to study the nature of the parameter space away from the domain wall, at a generic point of the 
two-dimensional [2] parameter space of V. As Figure [l] shows, the cusp splits the parameter space into two domains. 
Away from the domain wall, the three critical points of V are distinct. In the upper domain, two of the three roots 
are imaginary; there is only a single minimum. Below, they are all real. Notice that b —> —b has the same effect as a 
parity transformation, x — > —x. 

The physics here is stable. Near any value of a, there are two disconnected neighborhoods of b which support 
inflation. Perturbations of the scalar potential simply push these neighborhoods around. Inflation is possible as 
long as the physical parameter space includes the domain wall. This structure appears in the analysis of higher 
catastrophes. We shall see it embedded in the domain structure of A4 in section |IV| This implies that higher order 
corrections cannot exclude inflation unless they exclude a domain wall. This is the spirit of Catastrophe Theory. 

Apart from the neat description of parameter space, there are concrete, analytic results for inflationary physics. 
The separation of roots, (/3 — a), parametrizes many physical observables in a simple way, which we now investigate. 



B. A3 Inflation 



Very close to the inflection point, a, we expect to find cubic behavior of the potential. To that end, let y — x 
Substituting into Eq. (III.l), together with Eq. (III. 3 1 yields, 



a. 



V(y) 



1 



.or 



a)y 3 + X iy +(^ + c). 



(III.7) 



G 



The Two Domains of Cusp Parameter Space 




FIG. 1. A divide between two different regions of Cusp parameter space. 



Here, as in (II. 5 1, Ai represents a deviation in C from the domain wall. It is the slow roll parameter, e = Af/2. 
Eq. (II.7) gives number of e-foldings of inflation, 



N K = 



The Hubble parameter during inflation is related to the scale of inflation, which we may compute from (III. 7) 



Here V(8) is the final vacuum energy of the inflaton field. Finally, we come to the density perturbations. Integrating 
in Eq. ( II.6[ ) instead to xqobe where COBE perturbations leave the horizon gives -/Vtjobe, which is used to normalize 
the computed density perturbations to observations. Replacing x cn d and xcobe with N e and -/Vcobe, one finds that 
the density perturbations, 



A 2 



4tt 2 \dx/dN e J 



may be written as [8] 



A 2 - 



12tt 2 



W - <*) 2 N 4 e 



Immediately we sec that density perturbations scale parametrically with (/? — a), 



,2 V ( a ) AT A ta n2 



If, as in our universe, V(/3) ~ 0, we have the simple result, 



N 4 

A - = y °144^ a)6 - 



The physical observables discussed scale with (3 — a. As one may read from Eq. (III. 5 1, a and b vanish when a — (3 



The energy difference between inflation and the vacuum also vanishes. This is the consequence of V having a triply 
degenerate critical point at x — 0, that is, 



V = -x 4 + c. 
4 
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Notice that N e and A 2 ^ must be recalculated for this case. 

We close this section with Gilmore's [27] remark about the moniker "Catastrophe Theory". Had we taken the 
dual catastrophe, v4_ 3 with germ — x 4 /4, rather than A 3 , the same arguments would result in the function whose 
graph is reflected about the x axis. The three extrema phase would include a single, meta-stable minimum between 
two maxima. Upon passing through the anti-Cusp boundary, the system would lose metastability, tending towards a 
violent runaway. Catastrophe! 



IV. A A : THE SWALLOWTAIL CATASTROPHE 



Realistic models, like flux compactifications, typically involve more structure than an A3 potential affords. Indeed, 
the KKLT mechanism often includes a barrier to decompactification associated with SUSY breaking. This structure 
is only slightly more complex, and can be found in scalar potentials with A4 structure. In this section we describe 
inflation for A4 potentials along the lines of Section |III| We then apply these results to the study of Kahler moduli 
inflation in Type IIB String Theory. 



A. A4, Domain Structure 



The most general, local form for a potential with A4 behavior is 



V(x) 



A 



(IV.l) 



Here A may be positive or negative. Alternatively, one can take A > and connect to other models by a parity 
transformation, x — > —x. As with A3, the parameter space is divided into domains with different numbers of local 
minima. Again, we investigate inflation near doubly degenerate critical points. The same requirements of inflation 
and a local minimum forces the fourth critical point, 7, to be real, thus 



dV 

dx 



\(x- a) 2 {x- 0){x-i) 



Taking a derivative of Eq. (IV.l I and matching coefficients with Eq. (IV.2), we find 

2a + /3 + 7 = 0, 

and 



a/A = -(3a 2 + 2a/3 + /3 2 ), b/X = 2a(a + /3) 2 , c/A = -a 2 /3(2a + , 



(IV.2) 



(IV.3) 



(IV.4) 



Like A 3 potentials, the behavior of V depends on the sign of a. For a > parameter space has two domains: zero or 
one minimum. A domain with two minima opens only when a < 0. Demanding inflation and a metastable vacuum 
restricts are attention to the latter. Precise analysis of these regions can be accomplished with a trick from [27] . 
Under the scaling x — > C^j demanding V — > £ 5 , we find 



a -> C 2 a, b -> ( 3 b, 



C 4 c 



(IV.5) 



It suffices, then, to consider the cases a — 0, ±1 to understand the full behavior of the parameter space. The rest 
may be accessed via scaling. We shall consider the case of interest, a 
equations for the domain wall: 



T. Eliminating (3 in Eq. (IV.4), we have the 



b/X 



|a/A| 3 / 2 



-4a 3 + 2a, 



c/A 
WW 



3a 4 



(IV.6) 



a may be arbitrarily negative. To simplify (IV.6), one may simply redefine parameters to effectively take A = I. 
Fixing £, we have a cross section of the diagram, parametrized by the location of the degenerate root a. This is 
illustrated in Fig [2} 

At a generic point on this curve, V has an inflection point suitable for inflation. This is illustrated in Fig [3] The 
two corners correspond to triply degenerate critical points. One can easily check that the left-hand corner behaves 
like the A3 Catastrophe. The right-hand corner has A_3 behavior. At the top vertex corresponds to two, separated, 
doubly degenerate critical points: an Ai pair. The triangular, two minima region collapses when a = 0. Here, four 
extrema degenerate into a single critical point. At this point we find the quartically degenerate critical point, the A4 
germ. For a > 0, a parabola acts a domain wall between an unstable region and one with a single minima. We shall 
not consider this case further, but move on to physical considerations. 
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Various Domains of Swallowtail Parameter Space (a=-l) 




FIG. 2. A divide between two different regions of Swallowtail parameter space, for fixed a < 0. 
Swallowtail Catastrophe Model 

V(x) 




FIG. 3. The Swallowtail potential has a doubly degenerate critical point at generic points on the domain wall of parameter 
space. This is illustrated for the a = — I cross section of the three dimensional space. Note the similarity of structure to the 
typical Racetrack potential. 



B. A± Inflation 



We now investigate the properties of V near a domain wall. As with (III.7), let y = x — a, where a is again the 
inflection point. 



V(x) « V (^-(x - a) 3 + \i(x ~ a)) + V(a). 



(IV.7) 



The "scaling coefficient" of the x 3 behavior near the point of inflection is given by 



A 3 



-{a- P)(a-j). 
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An inflection point suitable for inflation exists when Ai = 0, which is precisely true if Eq. (IV. 2) holds. This 
expression can be simplified with Eq. (IV.3I. Applying this constraint is tantamo unt sh ifting the entire potential by 
a fixed (x). 



IV A 



This simplified the understanding of Swallowtail behavior in Section 
more transparent by forgoing this choice and leaving the location of the barrier, 7, explicit 
A3 is directly related to the density perturbations, 

A 2 _ V(a) 2 4 



From here on, the physics is 



12tt 2 



7V e 4 Z 2 (a-/?) 2 (a- 7 ) 2 



(IV. 



Here, N e denotes the number of e-foldings, and Z is a prefactor that depends on the canonical normalization of the 
inflaton. Again, the density perturbations scale with the separation of the roots of V'(x). 

Physically, V(a) corresponds to the scale of inflation and V(f3) is the cosmological constant. The energy scale of 
the barrier height, V(7) often carries other physical significance. For the example of Kahler moduli inflation discussed 
later, it is correlated to both the scale of SUSY breaking and decompactification of the internal space. Catastrophe 
Theory predicts scaling relations between them, as we now see. 

For our purposes, V(/3) is vanishingly small, as is Ai, for reasons mentioned earlier. From this, we can read off the 
scale of inflation 



V(a) = V { -^(f3-a)\ 



(0_a)*( 7 -a). 



(IV.9) 
(IV.10) 



The scale of the barrier height is, 



Vh) = V i -^[(P~a) 3 ~( 1 -aY 



(IV.ll) 



In short, the energy scales of the problem scale with both A3 and, more importantly, the separation of the roots. A 
large (f3 — a) corresponds to a large scale of inflation. With the interpretation of ^(7) as the scale of SUSY breaking, 
nearby roots imply that the scales of inflation and SUSY breaking are tied together. Thus, these scaling relations put 
the observations of Kallosh and Linde on precise, quantitative grounds. Additionally, we find that the scale of SUSY 
breaking can be separated from the inflation scale, and scales parametrically with the root separation (7 — a). 



V. RACETRACK INFLATION 



A. Parameter Space 

In this subsection, we analyze an example of swallowtail inflation: Kahler (volume) modulus inflation in type IIB 
string theory. We deal with a single Kahler modulus in the KKLT paradigm, keeping multiple Kahler moduli for 
future work. First we outline the parameter space of the physical input parameters that lead to successful inflation, 
with zero cosmological constant. Then, we comment on the scale of inflation and its relation to supersymmetry 
breaking. 

The main elements in the 4D Af — 1 effective theory are: (1) background flux induced Gukov-Vafa-Witten su- 
perpotential |25] . and (2) gaugino condensation on £>7-branes or Euclidean D3 instantons [32], which result in 
non-perturbative contributions to the superpotential. These two contributions to the superpotential are sufficient to 
stabilize complex structure moduli and the dilaton. Moreover, the Kahler moduli are stabilized in an AdS vacuum. 
An additional contribution to the scalar potential coming from anti-D3-branes then lifts the solution to a de Sitter 
vacuum 

The tree level Kahler Potential is 

K = -31og(T + f). 

In |12) the importance of a' and g s corrections to K for inflation was emphasized. These corrections can be 
parametrized schematically as, 
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In our analysis these additional parameters a% and £ are absorbed into the Swallowtail control space. They are 
examples of perturbations of V that do not change the domain structure of the Swallowtail Catastrophe. Thus, we 
no longer make reference to them explicitly. They may behave as nontrivial control parameters in multifield models. 
The effective superpotential is 

W = Wo + Ae- aT + Be- bT . (V.l) 

The imaginary part is stabilized by the nonperturbative breaking of the axionic shift symmetry. The flux contribution 
to the superpotential is given by Wo — f CY3 G^Afl. Gaugino condensation of pure Super- Yang-Mills on two stacks of 
D7 branes gives the non-perturbative superpotential contributions. The factors a and b are of the form 2tt/M, where 
M is the rank of the respective product in the gauge group. Also, A and B are coefficients in the non-perturbative 
terms whose precise dependence on moduli comes through one-loop threshold corrections to the D7 gauge kinetic 
function and higher curvature corrections on the world volume of the D7s. 
The TV = 1 supergravity scalar potential may be computed, 

V = e K (K IJ DjWDjW- 3\W\ 2 ) + (y ^ )2 , (V.2) 

where the final, "uplifting" , term breaks supersymmetry and is treated phcnomcnologically, as is usual in these 
constructions. Such a term may come from a -D3 brane trapped deep at the bottom of a warped throat |33l 135] - 

We now investigate a local range of the scalar potential which displays an inflection point and a local minimum. 
Clearly, the number of physical input parameters is large and treating all of them as free would entail studying a 
much higher order catastrophe. We therefore fix parameters that have little effect on the inflection point. Assuming 
that the gauge groups are fixed, we have four parameters at our disposal, Wo, A, B,C . While we will treat them as 
free, phenomenological parameters, one should not lose sight of their discrete origin in terms flux quanta. The goal, 
of course, is to use Catastrophe Theory to peel away model dependant details with the hope of determining what is 
natural about inflation in these models. 

The flux contribution \Wo\ 2 sets the scale of the potential. In KKLT constructions, low-scale supersymmetry 
requires \W \ be very small [T5]. This is usually done by assuming the magnitude of G^ 3 ' ) component of G 3 is 
suppressed relative to the primitive G^ 2,1 ) part by at least one part in 10 10 . We may remove Wo from the set of 
control parameters which influence the appearance of an inflection point, keeping it as a parameter that fixes the scale 
of the potential. This defines a new set of control parameters 

(Wo, A, B, C) -> (1, A, B, C) = (1, A * C_) . (V .3) 

Wo Wo w 

Supersymmetric AdS minima can occur, with C = 0, when 

x ~ i log \ A\ and 
a 

' log|4l- CV.4) 



a - 



b 



In the above, x is the real part of the Kahler modulus T. 

Before uplifting, the racetrack scalar potential behavior can be described by the Cusp Catastrophe. Of course, 
supersymmetry remains unbroken in the final vacuum. While nondegenerate critical points are possible here too, such 
a point would have negative vacuum energy. 

The uplifting term is positive definite, and will not contribute an additional minimum. With three parameters 
(A,B,C) and two minima, we are inevitably led to consider the Swallowtail Catastrophe, whose germ is x 5 . 

At this point, we can use the description of the parameter space of the Swallowtail Catastrophe, described in 



Section (IV A I, to understand the allowed regions. Here, we have to disting uish b etween the physical control param- 



eters A, B,C and the canonical control parameters a,b,c appearing in Eq. (IV. 1 1, in terms of which the swallowtail 
parameter space is cleanly described. Note that mathematically, there is no preferred set of control parameters, and 
the physical parameters are on the same footing as the canonical ones - the potential is a real function i3], has at 
most four local extrema and three parameters. It therefore must be described by the swallowtail catastrophe near 
the region relevant for inflation. 
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B. Scale of Inflation and Supersymmetry Breaking 



In this subsection, we continue our analysis of Racetrack inflation by investigating relations between the scale of 
inflation and the scale of supersymmetry breaking. 

The barrier height V('y) is roughly the scale of supersymmetry breaking [36 , 



V(7) ~ TO 3/2 



(V.5) 



V("f) vanishes if and only the uplifting coefficient C vanishes. This breaks the A4 structure down to A3. Eqs. (IV.9) 
and Eq. ( IV .111 suggest that low-scale supersymmetry entails low scale inflation. 

We are investigating the behavior of the Racetrack scalar potential, V in the vicinity of known extrema given in 
Eq. (TO). 



It is clear from Eq. (IV.9) that low-scale inflation is naturally obtained when the root separation is small. The 



scaling relations get more complicated with increasing complexity of the Catastrophe germ, but they are predictable 
and they scale paramctrically with the separation of the roots. However, from Eq. (IV. 8), it is clear that unacceptably 



small density perturbations are generic to this case unless Z, which varies parametrically with the Kahler modulus 
VEV at inflation, a is large. Z depends on the inverse Kahler metric; in the case of a single Kahler modulus, one 
obtains 



4a 2 

IT 



3/2 



(V.6) 



From Eq. (IV.8), we obtain 



(V.7) 



Typical values are A 2 ^ ~ 10 10 and N e ~ 50. For an intermediate scale of inflation, Vq ~ 10 28 . Then, for the correct 
value of density perturbations, we need 



a ~ O(10 2 - 10 3 ) . 

Note that the position of the inflection point is roughly given by its value prior to uplifting 

a ~ ~ log ™H • 



(V.8) 



(V.9) 



2 => Wo 

These large values naturally appear with magnetic flux threading the 4-cycles which the 



10 14 . Thus, Eq. (V.9) and Eq. (V.8) imply that for low-scale inflation, one requires 



We have V ~ \ W 
A ~ B ~ O(e 10 ° - e 

D7 branes wrap. The gauge kinetic function on the D7 branes is shifted by a piece /s that depends on the magnetic 
flux and the dilaton VEV 0. 

The struggle pointed out by Kallosh and Linde [36] was that the "uplifting" required by large scale inflation typically 
destabilize the potential. For nearby roots, one can obtain high scale inflation and low scale supersymmetry, without 
destabilization, by supplementing the model with additional constructions |31j . Another route is to attempt the 
separation within the LARGE volume scenario [16j 122 EH • This work obtained a hierarchy between inflation and 
SUSY breaking through a hierarc hy in f ield space between the modulus VEV where inflation occurs, and the position 
of the final vacuum[4 . . From Eq. ( IV. 1 1 ) it is clear that that (/3 — a) controls the scale of inflation and (7 — fi) controls 
the scale of SUSY breaking. Thus, any intermediate scenario is a priori possible. 

Three lessons come from the analysis of the inflection point. First, small root separation leads to similarity between 
the scale of inflation and SUSY breaking. This was seen in the example of volume modulus racetrack inflation. 
Second, if both SUSY breaking and inflation are near the intermediate scale 10~ 7 , then the inflection point needs to 
be located around a ~ O(10 2 — 10 3 ). Third, this implies that scale of inflation can be much larger than the scale of 
SUSY breaking for a -C j3 and j3 ~ 7. 



C. Racetrack Inflation: Numerical Results 



We now investigate the full scalar potential, Eq. (V.2) to demonstrate these ideas explicitly. After making our 
parametrization precise, we discuss the general general features and display the results in the plots below. 
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Low scale inflation requires large (x) and small root separations a — f3. One way to achieve this is involves turning 
on magnetic flux on the 4-cycles with the D7 branes wrap. The flux energy exponentially expands the internal volume, 
for details see [5]. 

We parametrize the magnetic flux by /s, and take 

A,B -> e afs A, e bfs B. 



It is easy to see that Eq. ( V.4) cause a shift x — >• x + fy. This construction allows us to absorb the entire effect of the 
parameter A into /j> 



Note that following the discussion after Eq. (V.9), it should be clear that /s is fixed by a and the scale of inflation 
given by Wo; specifically, the relation is 

fs ~ a + - log Wo . (V.10) 
a 

We may think of the family of potentials with the control parameters /, B, C as a bundle over the curve given by 
V. Because we are interested in the local behavior of such families of potentials, one should assign the swallowtail 
parameters, a, o, c as coordinates describing a section near a local trivialization of the bundle around some specific 
point. In other words, one should be careful to avoid large excursions in the space of Swallowtail control parameters. 

Concretely, we consider the fifth order Taylor polynomial as our Swallowtail function. To that end, we expand V 
around (x), 

5 

V(x) = J2 an{x - (x)) n + 0((x - (x)f). 
n=0 

The point (x^ should be fairly close to the critical point. A good starting point is to construct V around a known 



racetrack minimum, Eq. (V.4 1 



Restricting to the fifth Taylor polynomial (Jet space) is the simplest diffeomorphism consistent with the Swallowtail 
behavior, but has limited validity. A useful way to parametrize the faithfulness of this mapping comes from looking 
at 05 as a function of the Racetrack control parameters. Since the germ of V should be — we expect 05 < 0. One 
must be careful then to avoid portions of parameter space where 05 approaches zero[5]. 

A quick change of variable, 

y = (±^Ml _ ^1 (V 11) 

V (5o B )V8 5o 5 ' [VAl) 

gives the canonical form of the swallowtail potential, parity adjusted to agree with the Racetrack, 

1 r 1 o 1 9 

V = --y° - -ay s + -by 2 - cy + d. 



where 



50,50,3 — 2a| 
(5a 5 )8/5 ' 

25a 2 a5 — 15030405 + 4a| 
(5a 5 ) 12 / 5 ~ ' 
125aiag — 56020405 + 15030^05 — 3a| 



a 



(5o 5 ) 16 / 5 

25030405 + 125020405 — 625010405 



(5a 5 ; 



4 



We are then free to compare b and c to the domain structure of the Swallowtail parameter space. From a, b and c, 
one can use ( IV.4| and (IV. 3 1 compute the roots a,/3 and 7 



We investigated solutions near B = —0.0035, C — ^10 3 Wq and /s = 300. To fix the overall scale of the potential, 
we chose W = 2 x 10~ 10 . 

Fig|4]plots curves of various B parametrized by C with = 300. Greater B shifts the curves to the left. Inflation, 
which may occur near the triangular boundary, requires less uplifting parameter C . It generically corresponds to the 
decreasing separation of roots. Thus, B is bounded from above by the left vertex of the triangle: the Cusp point. 



13 



Here a = (3. B is bounded from below by the triangle's apex, which corresponds to the loss of a stabilized minima. 
In terms of roots, /? = 7. This behavior matches the full Racetrack potential. Smaller (more negative) values of B do 
not allow for metastable vacua. The same qualitative behavior occurs for any 

Fig |] plots curves of various /s parametrized by B with C = 2.7 x 10 22 . Increasing the flux, /s, pushes the 
curves leftward towards the Cusp point. Thus, the physical scale (greater than 30 TeV, say) of inflation bounds /s 
from above. The approximate fixed point for large B is approached in the region of a single minimum. As B is 
increased, the minima are stretched out and the barrier is flattened towards a decaying exponential behavior which is 
essentially independent of fe. Of course, SUSY is broken by the nonvanishing value of C, so a runaway is the limiting 
case here. Again, the Swallowtail parametrization faithfully reconstructs the qualitative and quantitative nature of 
the Racetrack potential. With such analytic control of the racetrack potential, we can easily generate solutions with 
zero cosmological constant which was a numerical challenge in earlier works. For low scale inflation, the cosmological 
constant needs to be vanishingly small so that it does not contribute to the soft supersymmetry breaking masses. 

Fig [6] shows an example of a curve of constant uplift parameter C, in the full, three dimensional Swallowtail 
parameter space. In the middle domain, the function possess two minima. Changing B changes the root separation. 
Smaller root separations mean less uplift is required to attain an inflection point. The curve pierces the domain wall 
at a value of B which gives an inflection point suitable for inflation. Further changes brings the function into a domain 
of a only a single minimum. Ultimately, the curve pierces the a second domain wall. Here both the minima are lost 
and the potential has runaway behavior for any x. 



Curves of Constant / £ in Swallowtail Control Space 

c 




FIG. 4. Curves of constant B parameter are plotted in the Swallowtail parameter space. Here the magnetic flux parameter is 
fixed at 300. As B is increased, the curves parametrized by C move to the right. Eventually they will pass the center of the 
triangle, where no stable vacua exist. As B decreases, eventually it will run into a cusp, requiring less and and less uplifting to 
reach the inflection point. 



D. Soft Masses and Control Parameters 

When the inflaton also participates in supersymmetry breaking in the vacuum of the theory, inflationary parameters 
directly influence low-energy particle physics, such as the supersymmetric spectrum. In this subsection, we outline 
the dependence of soft masses on the position of roots in the inflaton space, as well as the root separation. 

Firstly, we recap the M = 1 effective theory. The visible sector is assumed to be sequestered from the supersymmetry 




FIG. 5. Curves of constant magnetic flux are plotted in the Swallowtail parameter space. The value of B where each curve 
intersects the triangular region, and inflection point suitable for inflation exists on the Racetrack Potential. As magnetic flux 
increases, the curves come closer to the Cusp point, implying a lower ratio of inflationary and cosmological constant scales. 
This is all for C = 2.7 x 10~ 22 . 



breaking brane, [33]. The model data are, 

k = -3io g (T + f) + H T , r)** **, 

i 

W = W + Ae~ a ^ + Be- b ^ + ( '. A, , ; <I»,<I» ,<1>,.. 

fg = T-fr, (V.12) 

where $i are the visible sector superfields and Zi(T,T) = 1/(T + T) ni . rii is the modular weight of the superfields - 
either zero or one if they live on D7 or D3 branes, respectively, /s is the contribution to the gauge kinetic function 
coming from magnetic flux on the D7 branes. Recall that our effective variable is x = Re(T). 

Since the vacuum has zero cosmological constant, we have only moduli and anomaly mediated contribution to the 
soft SUSY breaking masses at the GUT scale. For more details of mirage mediated SUSY breaking, see [TM2T1 12rj] . 
The modulus contributions to gaugino masses, trilinear couplings and sfermion masses are given by 

AI Q = F T d T logRe(f g ), 
m\ = -F T F T d T d T log(e- K "^Z l ), 
Aijk = F T d T log(e- Ko Z i Z j Z k ) . 

(V.13) 

Where Kq is the Kahlcr potential without contributions from the superfields, $i. As we have discussed previously, 
the SUSY breaking scale is similar to the barrier height V(7), which in turn is similar to the inflationary scale V{a) 
for small root separation. Thus, 

F 2 - U( 7 ) - V(a) 

for a ~ P ~ 7 . (V.14) 



It follows from Eq. (IV.8) that 
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Curve of Constant B parameter in Full Swallowtail Parameter Space 

-2.0 





FIG. 6. Curve of constant C (thick solid line) in the three dimensional Swallowtail parameter space. As the parameter B is 
changed, (j3 — a) also changes. When the curve pierces the inner surface, B is suitable for inflation. As the roots approach 
each other the uplifting is too strong and the function is in the domain with only a single minimum. For a large change in B, 
the curve passes through the second domain wall where V has runaway behavior for any x = Re[T]. 



It is instructive to see what this means for the modulus contribution to the GUT scale gaugino mass, for example. 
Using Eq. (V.15I and the expression for the gauge kinetic function, one obtains 

u = 8ir\A n \a 3 , , 

3iV2|(/3-a)(7-«)|(/3-/ S ) ' 1 ' ' 

Note that we have used x = (3 as the location where the gauge kinetic function is calculated, giving f g = /3 — fy. 
The anomaly contribution to the gaugino mass at the GUT scale is given by 

Manomaly = > (V-17) 

where b is the beta function coefficient for the relevant gaugino, and g is the gauge coupling. The gravitino mass and 
the modulus contribution Mq are related by a parameter a m i r that controls the relative strengths of modulus and 
anomaly contributions 



M log(M p /m 3/2 ) 

Finally, a m ir can be obtained as a function of the position of the inflaton [5] 

amir _ 32 y p 



(V.18) 



■ (v - 19) 
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Putting everything together, we obtain that the total GUT scale gaugino mass is given by 



M 



8ir\A n \a 3 



total 



3jV e 2 |C8-a)( 7 -a)|03-/ s ) 



1 



167T' 



log(M p /m 3/2 ) 



(V.20) 



10" 



which follows from the value of Ar - 10" 5 and Eq. This 



First, we can fix the scale Mq 
leads to gaugino masses in the 100 — 1000 GeV range after RG running. The precise value would depend on the exact 
separation of roots. Second, the relative importance of modulus versus anomaly mediation can be controlled by the 
root position j3, as is clear from Eq. (V.19I. Note that as mentioned before, is not an independent parameter but 
is given by Eq. (V.10). 



This suggests that if we can identify a m i r the ratio of the mirage to anomaly mediation contribution and Mq at 
the LHC [25], we can determine two differences among the three roots, a, (3 and 7. Cosmology data may determine 
the scale of inflation, and therefore combining the data from the LHC and cosmology, we may be able to identify the 
nature of the potential responsible for inflation and SUSY breaking. 



VI. CONCLUSION 



Many attempts to incorporate inflation within the framework of particle physics and string theory involve inflection 
points. While many models exist with markedly different motivations, we have taken the point of view that they 
are all related to universal archetypes given by Singularity Theory, and tried to extract physical information in a 
model-independent way. Our approach gives 

(i) A neat classification of the kinds of local potentials that one can obtain for a given number of physical input 
parameters affecting inflation. 

(ii) A full description of the parameter space where one can expect inflationary solutions. The space of couplings or 
control parameters is divided into domains with different qualitative behaviors of the scalar potential. Near a domain 
wall, a function possesses an inflection point suitable for inflation. Perturbations away from the degenerate critical 
point determine the number of e-foldings of inflation. These domain walls depend on the separation of the roots of 
V'(x) in a precise way. These features are stable against perturbations as long as a domain wall is not crossed. 

As mentioned in the text, the number of e-foldings N e parametrizes the unfolding of the singularities associated 
with degenerate critical points. To get a physically reasonable number e-foldings, the parameters of the potential 
must be very close to a domain wall. Universal behavior in Warped D-Brane inflation was characterized in terms 
of e-foldings in [BJ. It is possible to generalize these results to all models using inflection points, and we have found 
preliminary evidence in agreement with [6], which will be presented in a later publication. 

(in) The fact that domain walls depend on the separation of the roots of V'(x) in a precise way turns into scaling 
relations for physical observables such as density perturbations and the energy scales of inflation and SUSY breaking. 
For low-scale inflation with metastable vacua, this analysis quantitatively explains Kallosh and Linde's bound between 
the inflationary and SUSY breaking scales. Low-scale SUSY breaking with high-scale inflation may be achieved with 
large root separation, regardless of the origin of the scalar potential. The scaling relations even appear in computations 
of the soft SUSY breaking masses. Successful measurement of SUSY masses at colliders will allow correlations between 
particle physics data and cosmology data through archetype models such as the Swallowtail. It could also rule out 
inflection point inflation altogether. 

In short, we have taken inflection point inflation (like phase transitions) as an application of Catastrophe Theory. 
The application gives both qualitative and quantitative information. 

A singularity theoretic approach to inflation may also be useful for understanding the behavior of potentials as 
the VEV of the inflaton approaches the Planck Scale. In this limit, new, higher dimensional operators increase the 
dimension of the space of control parameters. Since lower order catastrophes nest inside higher order ones, inflation 
is still viable, but the control parameter space must be appropriately analyzed. It may be possible, then, to organize 
the Planck suppressed operators and make contact with particle physics. [7] is a first step in this direction. 

As previously mentioned, N e parametrizes the unfolding of the singularities associated with degenerate critical 
points. To get a physically reasonable number e-foldings, the potential's parameters must be very close to a domain 
wall. Universal behavior in Warped D-Brane inflation was characterized in terms of e-foldings in [B] . It is possible to 
generalize these results to all models using inflection points, and we have found preliminary evidence in agreement 
with [BJ, which will be presented in a later publication. 

This work has focused on ^4-type singularities. Multi-field inflation brings about the possibility of singularities of 
D- or E-type. Indeed, searches for multi-field models of Kahlcr Moduli inflation are numerically challenging. Yet 
such models seem to naturally incorporate B-type singularities. Similar techniques to those in this paper are readily 
available, and may have interesting consequences for the study of non-gaussianities in the power spectrum. We shall 
report on these details in future work. 



17 



VII. ACKNOWLEDGEMENTS 



This work is supported in parts by the DOE grant DE-FG02-95ER40917. We thank R. Allahverdi for useful 
discussions and the participants of the Aspects of Inflation workshop at the Mitchell Institute for stimulating questions. 



Appendix A: Critical Points and Symmetry Breaking 



1. Symmetries of Critical Points 



Catastrophes of a single variable stem from the germ, x n /n. Under Arnold's classification 1 J , it is an A n ^i critical 
point. The fully perturbed germ is then 

1 n— 2 

f n (x) = -x n + Va m x m . (A.l) 
n ' 

One way to understand the origin of the A n -\ symmetry of the critical point is to interpret /„ as the characteristic 
polynomial of some n x n matrix, M, in the adjoint representation of A n _\ (that is, SU(n)). The control parameters 
a,i are then specific sums of products of traces of products of M, as one might be familiar with from studying the 
Cayley-Hamilton Theorem. For example, 

a n -2 = i[(trAf) 2 -tr(M 2 )]. 

Other important examples arc ao = dctM, and a„_i = — trM = 0. The latter is a generic property of "Gell-Mann" 
matrices of A n -\. Notice that these coefficients, by the cyclic property of trace, Oj are manifestly invariant under a 
A n -i group action. 

M -> g*Mg, 



for A n _i element g. Generally, M is a vector in the complexified Lie Algebra of A n -\. The result of this analysis is 
that the control parameters of Eq. (A.l) are maps from this linear space to C. Typically we restrict to cti <E R, but 
some nontrivial phases may be important in certain models of inflation, and as such we get the added bonus of being 
able to understand them. In either case, the complexification is still important, for if 



with Ai being the Gell-Mann matrix associated with a basis element of the Cartan subalgebra of A„_i, then 



n-1 

a n -2 = ~J2 c l 



That Cj € C is important to ensure the possibility that a„_2 can be either positive or negative. Aside from filling the 
entire Catastrophe parameter space, a„_2 = is generically located on a domain wall. The upshot all this is that the 
domain structure of control parameter space maps nicely into the Weyl chambers of the weight space of A n _i. 



2. Breaking Symmetries of Critical Points 

Breaking a Catastrophe of type A n _\ down to type A n _2 is a simple matter of analyzing germs, 

n x 



n n n — 1 

Here t parametrizes the breaking. Demanding that t ^ is equivalent to demanding that the matrix M have a 
nonvanshing trace, breaking the A n -\ symmetry. We then find ourselves with a germ of the A„_2 catastrophe. 
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For the polynomials we have been considering, this merely amounts to shifting x — > x + t. In the case of the Cusp 
Catastrophe, this amounts to moving along the curve depicted in Fig[l] For a general n, this breaking can be extended 

— -> — + — x"- 1 + -^x n ~ 2 + ... + -x 3 . 
n n n — 1 n — 2 3 

Any further breaking will completely morsify the function to which this Catastrophe belongs. Tuning the parameters 
away from zero effectively "blows up" the various spheres in parameter space, circles in this case, associated with the 
mapping singularity. 

The careful reader may complain that we have simply renamed the control parameters. For the case of a polynomial 
functions, this is true, modulo the shift. It is in this sense that polynomials are the archetypes for more complicated 
functions, and our philosophy is that inflation need only be understood in terms of said archetypes. The point is that 
this "Catastrophic" behavior exists so long as the function is smooth in the relevant domain. Put another way, a 
simple Taylor expansion is not always the appropriate diffcomorphism to analyze this behavior. 

Concretely, let V(x) be a function whose parameter space allows for a quartically degenerate critical point — 



Swallowtail behavior. As in Eq. (II. 3), we have 

^ = (x - a){x - (3)(x - 7 )(z - d)f(x), (A.2) 

where f(x) is nonvanshing for any x. If each root is different, we have a morse function with four local extrema. If 
a = 5, say, we'd have a doubly degenerate critical point. A typical flux compactification scalar potential may have 
f(x) ~ e~ x /x n , which means that V is certainly not a polynomial! 

If the Swallowtail parameters are quite small, a fifth order Taylor polynomial will describe the functional behavior 
quite well. As the roots begin to separate, the Taylor series will loose fidelity. Despite this, Eq. |A.2| still guarantees 
the structural behavior. The domain walls in parameter space still exist, but the mapping from control parameter 
space to the simple Swallowtail model becomes complicated, perhaps intractably so. 

Physically, however, this shouldn't concern us. What matters is simply that a doubly degenerate critical point, a 
local minimum and a barrier to runaway exist. From this perspective, we can simply look for the catastrophe 
relevant for inflation. Let y = x — a. The fifth order Taylor polynomial becomes 

V(x + a) « Voi^y 5 + \y A + S -y 3 + ^y 2 + My) + V(a). 

The parameters t and s now contain the data about the barrier and metastable vacuum. If we choose a to be the 
inflection point, then, 

V(x + a) = V (^y 3 + An/) + V(a) + 0(y 4 ), 

which is is all the data we need to investigate the number of e-foldings, as we now demonstrate. Ai serves as the slow 
roll parameter e. Because of the small field excursion, the number of e-foldings will be given very accurately by 



dy 
sy 2 + Xi 



The rest of the physical properties follow from similar considerations. 

In all, a singularity theoretic approach to studying inflation is particularly well suited to embedding it within a 
larger physical theory. Existence of inflation depends on existence of domain walls in allows regions of parameter 
space. Relationships between inflation and other phenomena, like SUSY breaking are illuminated. This behavior is 
necessarily universal. 



[1] A simple example of a nontrivial restriction of the parameter space would be two widely separated degenerate critical 
points. While they do descend from the unfolding of a critical point of higher degeneracy, it may be physically simpler 
to consider them as isolated Catastrophes with their own parameter spaces with some appropriate matching conditions. 
This is reasonable because Catastrophe theory involves the local behavior of functions. Indeed, the intuition gained from 
"gluing" is helpful when considering more complex Catastrophes, such as those that involve multiple variables. 
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[2] This is certainly true in the slow roll approximation, Eq. ( II. 6 1. Interestingly, this remains true when a broader set of initial 
conditions is considered. This is be explored in future work. 

It was shown in [T5] how to set all the parameters of these sorts of constructions to real numbers. 

The basic idea is that if some amount of background sources like radiation is present, an attractor solution to the evolution 
equations can guide the field to the global minimum at x = p without overshooting, even when the extrema are widely 
separated. 

Although this can be often remedied by a simple diffeomorphism: shift x. A more complete treatment is discussed in 
Appendix |A"| 
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